Exclusion of Tiny Interstellar Dust Grains From the 

Heliosphere 

J. D. Slavin*, P. C. Frisc^, J. Heerikhuisen**, N. V. Pogorelov**, H.-R. Mueller*, 
W. T. Reach 1 , G. P. Zank**, B. Dasgupta** and K. Avinash** 

""Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 (USA) 
^Dept. Astronomy & Astrophysics, U. Chicago, 5640 S. Ellis Ave., Chicago, IL 60637 (USA) 
**U. Alabama-Huntsville, Center for Space Plasma and Aeronomic Research, Huntsville, AL 35899 
^Dept. of Physics and Astronomy, 6127 Wilder Lab , Dartmouth College, Hanover, NH 03755 
§ Infrared Processing and Analysis Center, Caltech, Pasadena, CA 91125 

Abstract. The distribution of interstellar dust grains (ISDG) observed in the Solar System depends on the nature of the 
interstellar medium- solar wind interaction. The charge of the grains couples them to the interstellar magnetic field (ISMF) 
resulting in some fraction of grains being excluded from the heliosphere while grains on the larger end of the size distribution, 
with gyroradii comparable to the size of the heliosphere, penetrate the termination shock. This results in a skewing the size 
distribution detected in the Solar System. 

We present new calculations of grain trajectories and the resultant grain density distribution for small ISDGs propagating 
through the heliosphere. We make use of detailed heliosphere model results, using three-dimensional (3-D) magnetohydrody- 
namic/kinetic models designed to match data on the shape of the termination shock and the relative deflection of interstellar 
H° and He flowing into the heliosphere. We find that the necessary inclination of the ISMF relative to the inflow direction 
results in an asymmetry in the distribution of the larger grains (0.1 jUm) that penetrate the heliopause. Smaller grains (0.01 
jim) are completely excluded from the Solar System at the heliopause. 
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INTRODUCTION 

Up to one percent of the mass of the interstellar cloud 
surrounding the heliosphere is carried by dust grains that 
interact with the heliosphere [1, 2]. The flow of inter- 
stellar material past the Sun at 26.4 km s _1 drives large 
interstellar dust grains into the heliosphere, while small 
grains are diverted around the heliosphere by the inter- 
stellar magnetic field (ISMF). The heliospheric trajecto- 
ries of intermediate sized grains can be complicated and 
depend on the solar wind magnetic field and thus solar 
cycle phase. Observations of ISDGs in the solar system 
by the Ulysses, Galileo and Cassini spacecraft show that 
the density of smallest grains are deficient in the inner 
heliosphere when compared to the nominal Mathis et al. 
"MRN" power law size distribution [e.g. 3, 4]. For in- 
stance, the density of grains of mass 10 -14 g (a ~ 0.1 
jUm) is reduced in the inner heliosphere by a factor of 
~ 90 below the MRN predictions and 10" 15 g (a ~ 0.05 
jUm) grains are deficient by three orders of magnitude. 
At the same time large (a ~ 1 /im) grains are absent 
from the MRN distribution, but are abundant in the in- 
flowing ISDGs. While several alternative grain size dis- 
tribution models exist [e.g., 5, 6, 7], the constraints of 
staying within the limits of cosmic abundances and ex- 
plaining the interstellar extinction curve lead in all cases 



to the presence in the models of small grains that are ab- 
sent from the observed distribution and a cutoff on the 
large grain size end that is below that of the observed 
distribution. The deficiencies of small grains may result 
from small grains being deflected around the heliosphere 
because their high charge-to-mass ratios cause them to 
couple tightly to the ISMF. On the other hand, the cir- 
cumheliospheric interstellar medium (CHISM) may sim- 
ply be deficient in small grains (balancing the overabun- 
dance of large grains). To evaluate which is the case re- 
quires calculations of the grain distribution as a function 
of grain size in the context of realistic heliosphere mod- 
els. With recent data on the shape of the heliosphere from 
the crossing of the termination shock (TS) by Voyagers 
1 and 2 [8], and on the deflection of inflowing interstel- 
lar H° relative to He [9], such models face tighter con- 
straints than ever before. The inferred asymmetry of the 
heliosphere can be explained by the orientation of the 
ISMF relative to the direction of the cloud- Sun relative 
motion. The ISMF required to fit the heliosphere data, in 
turn, has direct implications for ISDG trajectories. 

In this paper we present new calculations of grain tra- 
jectories and density distributions for small grains that 
are primarily excluded from the inner heliosphere. Small 
grains experience enhanced charging rates due to the 
ejection of secondary electrons in the hot 10 5 — 10 6 K 



plasma between the termination shock and heliopause 
[10]. We calculate grain trajectories in 3-D, based on 
a self-consistent steady state MHD-kinetic heliosphere 
model with grain charge calculated at each location along 
the path according to the plasma conditions and radiation 
field. (We note that the recent Voyager 2 results [11] on 
the lower than expected plasma temperature in the inner 
helio sheath are not incorporated in the heliosphere model 
we use. We expect that this could result in somewhat 
lower dust charging in this region than we calculated, 
though the detailed effects of charging by thermal par- 
ticles vs. PUIs has yet to be fully explored.) The model 
of plasma density and magnetic field in the heliosphere 
at each point of space properly accounts for the asymme- 
tries that result from the angle between the ISMF and in- 
terstellar gas/dust inflow direction, and ISMF and eclip- 
tic plane [12, 13]. When combined with recent models of 
grain charging [14] and a realistic ultraviolet (UV) radi- 
ation field, the Lorentz force at each point in space can 
be used to calculate the grain trajectory. Total 3-D grain 
densities are then calculated for a sample of interstellar 
grain sizes based on the total gas-to-dust mass ratio in 
the CHISM and ISDG models. Many previous calcula- 
tions have been made of grain interactions with the he- 
liosphere [e.g. 15, 3, 10, 2] but ours is the first to use a 
realistic heliosphere model, including distortion due to 
an ISMF orientation that is consistent with recent obser- 
vational data. 



GRAIN TRAJECTORY CALCULATIONS 
Heliosphere Model 

The grain trajectory calculations require a 3-D model 
that self-consistently includes the interaction and charge 
exchange between interstellar H and the solar wind 
plasma, and that is based on realistic boundary condi- 
tions for the heliosphere as set by the surrounding in- 
terstellar cloud. The 10 AU difference in the termina- 
tion shock distances seen by Voyager 1 versus Voyager 2 
suggests an asymmetry in the outer heliosphere topology 
that can be explained by the tilt of the ISMF with respect 
to the upwind direction. We therefore select a heliosphere 
model that reproduces the TS asymmetry, the hydrogen- 
helium offset of ~ 5°, and the angle of the ISMF with 
respect to the ecliptic plane as indicated by the H-He off- 
set and data on the ISMF near the Sun. We use the he- 
liosphere of Pogorelov et al. [12, 13] and Heerikhuisen 
et al. [16], which is based on coupled 3-D MHD (ion) 
fluid and kinetic (neutral) particle code including the ef- 
fects of charge exchange to calculate the plasma and neu- 
tral particle properties. The model parameters of the solar 
wind and inflowing interstellar material are given in Ta- 




FIGURE 1. 3-D MHD/kinetic heliosphere model used in 
these calculations [12, 13], showing the heliopause boundary, 
magnetic field lines and the Voyager 1 and 2 trajectories. The 
viewpoint is from interstellar space looking towards the down- 
wind. The heliopause is colored according to the strength of the 
magnetic field. The bending of the magnetic field lines around 
the heliopause dominates the trajectories of tiny dust grains in- 
teracting with the heliosphere, and is responsible for the dust 
trajectories and hence the enhanced dust densities. 



ble 1 (where n(H°) is the neutral H density and ft(p + ) is 
the proton density). The solar wind in the model carries a 
Parker spiral type magnetic field and has v(H + )sw = 450 
km s" 1 . The ISMF strength is Bi SM = 3 jllG with the 
direction 30° away from the ISM inflow direction and 
pointed toward the southern hemisphere, and lying in the 
H deflection plane for this model [13]. The heliosphere 
model uses a coordinate system with the x-axis perpen- 
dicular to the ecliptic plane, the z-axis in the ecliptic 
plane towards upwind and perpendicular to the x-axis, 
and the y-axis in the ecliptic plane to form a righthanded 
coordinate system. The Sun is located at x = y = z = 0. 
Fig. 1 shows the ISMF lines bending over the heliopause. 



Grains and Grain Charging 

We illustrate the deflection of small grains around the 
heliosphere by calculating the charge and trajectories for 
spherical grains of radii 0.1 /im, and 0.01 /im, composed 
of astronomical silicates. The assumption that the small 
grains are silicates is supported by observations and anal- 
ysis indicating that the gas phase carbon abundance in the 
CHISM is high [supersolar in fact 17] leaving none to be 
locked up in grains. Grain charging is calculated using 
the model and grain-charging code of Weingartner et al. 
[14], combined with plasma densities and temperatures 
from the heliosphere model, and the modeled extreme- 



TABLE 1. Parameters of Heliosphere Model with Positive Magnetic Polarity at the NEP 





Velocity 


«(p+) 


»(H°) 


Temperature 


Magnetic Field 




km s _1 


cm J 


cm J 


K 




ISM at "oo" 


26.4 


0.06 


0.15 


6,527 


3.0 


Solar Wind at 10 AU 


454 


0.07 





2,750-14,280 


0.4-3.3 




FIGURE 2. Example of spatial trajectories for 0.1 jUm grains (left) and 0.01 jUm grains (right). The upwind direction is toward 
the back plane while the ecliptic is parallel to the bottom plane. Note the convergence of trajectories in the top and bottom, with 
most 0.01 fim grains deflected to low latitudes. The 0.01 jirni grains are entirely deflected around the heliosphere. 



UV/far-UV photon fluxes as discussed below. 

The primary mechanisms affecting the charge of 
grains near and in the heliosphere are the photoelectric 
emission caused by absorption of photons from the UV 
radiation field, and the sticking of electrons and ions that 
collide with the grains. Additional grain charging occurs 
in the hot subsonic solar wind plasma in the inner he- 
liosheath between the solar wind termination shock and 
the heliopause because of secondary electron emission 
produced by collisions of high energy electrons with the 
grains. 

The radiation field used in our calculations is com- 
prised of the constant interstellar radiation field from 
Gondhalekar et al. [18], and the Solar flux data taken 
from the TIMED/SEE project (see http : //lasp . 
Colorado . edu/see/see_data . html). The 
TIMED/SEE project creates spectra using a combina- 
tion of measured and modeled fluxes. The photon flux 
long ward of 13.6 eV is dominated by the interstellar 
radiation field, while the extreme ultraviolet radiation 
field is dominated by the Sun. The interstellar and solar 
fluxes have equal effect in charging grains at « 100 AU. 



Grain Trajectories 

The grain trajectory calculations begin with the grains 
upstream of the heliosphere and traveling with the inter- 
stellar gas. However, the grains are likely to have been 
accelerated by turbulence in the interstellar magnetic 



field and thus should have some initial velocity relative 
to the gas. Taking guidance from Yan et al. [19], we give 
the grains an initial 3 km s _1 speed (roughly the gas tur- 
bulent velocity observed for the Local Interstellar Cloud) 
relative to the gas, and with a net direction that is perpen- 
dicular to the magnetic field but is otherwise randomly 
directed. 

We integrate the equations of motion for the grains, 
calculating at each point the charge on the grain depend- 
ing on its position (and properties of the surrounding 
plasma) and including the Lorentz forces, radiation pres- 
sure and gravity of the Sun. For the density distribution 
we take over 10 6 grain trajectories for grains starting in a 
grid in a plane far (~ 900 AU) from the Sun in the ISM. 

The forces acting on the grains include gravity, radia- 
tion pressure, and the Lorentz force from the motion of 
charged grains relative to the magnetic field lines. The 
equation of motion is given by 

5 g + (l- J S)^x g -^((f g -v gas )x J B)=0 (1) 

where 3c g , q, and m are the velocity, charge and mass of 
the grain, j3 is the ratio of radiation pressure to gravita- 
tional attraction, B is magnetic field, and v gas is gas (solar 
wind or ISM) velocity. Although radiation pressure is not 
significant for spherical grains in the outer heliosheath 
that obey Mie scattering, Kimura and Mann [20] have 
shown that radiation pressure acting on fluffy grains re- 
sults in non-radial torques that may significantly perturb 
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FIGURE 3. Dust density distributions (cm -3 ) in a slice in the z = plane (perpendicular to the IS flow direction) for radius 
0.1 jim (left figure) and 0.01 jim (right figure) silicate grains. (See text for coordinate system.) These views are centered near the 
upwind direction, and ecliptic longitude increases towards the left. The inner nearly circular contour is the solar wind termination 
shock, while the outer contour is the heliopause. The a ~ 0. 1 jtim grains are partially diverted at the heliopause, but some grains 
penetrate inside the TS. All 0.01 jirni grains, which have small gyroradii due to a high q/m ratio, are excluded from the heliosphere. 
The colors indicate the dust density relative to its value in the pristine ISM. Note that the plotted density range is much larger for 
the 0.1 jirni grains. 
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FIGURE 4. Spatial densities of grains (cm -3 ) for slices in the xz plane. The right figure (0.01 jUm grains) shows the full exclusion 
of tiny grains from the heliosphere. The left figure (0. 1 jirni grains) shows the diversion of slightly larger grains able to partially 
penetrate the inner heliosheath. The concentration of grains near the North Ecliptic Pole results from the deflection of particles in 
the inner heliosheath where they are sensitive to the magnetic polarity of the solar wind. 



the grain motion in the outer heliosphere. This effect is 
not included here. 

We show sample grain trajectories for silicate grains of 
radii 0.01 jim and 0. 1 jim in Fig. 2. The resulting density 
distributions of grains (cm -3 ) are shown for a slice in the 
plane at z — that is perpendicular to the interstellar flow 
direction (Fig. 3), and for a meridian slice through the 
heliosphere nose (Fig. 4). 

These simulations of grain locations in the extended 



3-D volume were performed with the University of Al- 
abama, Huntsville, CSPAR cluster. 

RESULTS 

This study shows that the distortion of the local interstel- 
lar magnetic field by the heliosphere causes large-scale 
asymmetries in the spatial distribution of tiny interstellar 



grains flowing past the Sun and within 1000 AU. Grains 
of radius ~ 0.01 jim appear to be completely excluded 
from the heliosphere. The trajectories of these grains will 
be sensitive to the direction of the local interstellar mag- 
netic field, and relatively insensitive to the solar 22-year 
magnetic activity cycle since the grains do not penetrate 
the heliopause. In Fig. 5 we show the column densities 
of these grains within 400 AU of the Sun, and outside of 
the heliopause. The ISMF distortion by the heliosphere 
is imprinted on the distribution of these small grains. 

Grains as small as 0.1 jim can partially penetrate the 
termination shock. The trajectories of 0.1 jim grains 
are initially influenced by the interstellar magnetic field. 
Once inside of the heliosphere, the solar wind magnetic 
field and solar cycle become significant [e.g. 2]. The 
high charge-to-mass ratio for these grains is enhanced 
by secondary electron ejection in the inner heliosheath, 
increasing the coupling between 0.1 jim grains and the 
magnetic field in this region. 




FIGURE 5. Column densities (cm -2 ) of 0.01 jim grains 
deflected around the heliopause, out to a distance of 400 AU 
from the Sun. The values on the color bar show the true column 
density (cm -2 ) for the assumptions that the gas-to-dust ratio 
is 100 and all the grains are 0.01 jim. The upwind direction 
is at the center, with 0° corresponding to ecliptic coordinate 
£-255°. 
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